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■ Abstract 

We give an introduction to the heat kernel technique and ^-function. 
Two applications are considered. First we derive the high temperature 
asymptotics of the free energy for boson fields in terms of the heat kernel 
expansion and ^-function. Another application is chiral anomaly for local 
iy*~) ■ (MIT bag) boundary conditions. 
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In this paper we give an introduction to the technique of the heat-kernel expan- 
sion and ^-function regularization. The heat kernel became a standard tool in 
calculations of the vacuum polarization, the Casimir effect and quantum anoma- 



lies. In the case of non-trivial background fields, and especially in the presence 
of boundaries or singularities, the heat kernel technique seems to be the most 
adequate one for the analysis of the one- loop effects (see [I] for a recent review). 

A possible application of heat kernel methods to the evaluation of high temper- 
ature asymptotics of the free energy in the presense of boundaries or singularities 
has not received due attention yet. At high temperatures the expansion of the 
free energy can be determined by the heat kernel expansion and ^-function 

Chiral anomaly, which was discovered more than 30 years ago [3], still plays 
an important role in physics. On smooth manifolds without boundaries many 
successful approaches to the anomalies exist The heat kernel approach to 
the anomalies is essentially equivalent to the Fujikawa approach |Hj and to the 
calculations based on the finite- mode regularization jOj, but it can be more easily 
extended to complicated geometries. The chiral anomaly in the case of non- 
trivial boundary conditions (MIT bag boundary conditions) has been calculated 
only recently [Zj. 
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The paper is organized as follows. In Sec. 2 we give an introduction to the 
formalism of the heat kernel and heat kernel expansion. Also we introduce a zeta 
function and calculate the one-loop effective action in terms of the zeta function. 
In Sec. 3 we consider two examples. First we derive the high temperature expan- 
sion of the free energy for boson fields in terms of the heat kernel expansion and 
^-function. Then we discuss a chiral anomaly in four dimensions for an euclidean 
version of the MIT bag boundary conditions [SJ. 

2 Spectral functions 
2.1 Heat kernel 

Consider a second order elliptic partial differential operator L of Laplace type 
on an n-dimensional Riemannian manifold. Any operator of this type can be 
expanded locally as 

L = -(g^d^ + a°d a + b), (1) 

where a and b are some matrix valued functions and g^ v is the inverse metric 
tensor on the manifold. For a flat space g^ v = 5^. 
The heat kernel can be defined as follows: 

K(t;x;y;L) = (x\ exp(-tL)\y) = ^ <f>{(x)4> x (v) exp(-*A), (2) 

A 

where cj)\ is an eigenfunction of the operator L with the eigenvalue A. 
It satisfies the heat equation 

(d t + L x )K(t;x;y;L) = (3) 

with an initial condition 

K{0;x;y;L) = 6{x,y). (4) 

If we consider the fields in a finite volume then it is necessary to specify 
boundary conditions. Different choices are possible. In section 3.1 we will consider 
the case of periodic boundary conditions on imaginary time coordinate, which are 
specific for boson fields. In section 3.2 we will study bag boundary conditions 
imposed on fermion fields. If the normal to the boundary component of the 
fermion current ^ Tn^ vanishes at the boundary, one can impose bag boundary 
conditions, a particular case of mixed boundary conditions. We assume given two 
complementary projectors Il±, II„ + Il + = / acting on a multi component field 
(the eigenfunction of the operator L) at each point of the boundary and define 
mixed boundary conditions by the relations 

n_vu = o, (v n + s)n + v>U = o, (5) 
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where S is a matrix valued function on the boundary. In other words, the compo- 
nents satisfy Dirichlet boundary conditions, and H+i[> satisfy Robin (modi- 
fied Neumann) ones. 

It is convenient to define 

x = n + - n_ . (6) 

Let {ej}, j = 1, . . . , n be a local orthonormal frame for the tangent space to 
the manifold and let on the boundary e n be an inward pointing normal vector. 
The extrinsic curvature is defined by the equation 

= r™ fe , (7) 

where T is the Christoffel symbol. For example, on the unit sphere S*™ -1 which 
bounds the unit ball in R n the extrinsic curvature is L ab = 5 a6 . 

Curved space offers no complications in our approach compared to the flat 
case. Let R^ upa be the Riemann tensor, and let R^ v = R° '^ va be the Ricci tensor. 
With our sign convention the scalar curvature R = R£ is +2 on the unit sphere 
S 2 . In flat space the Riemann and Ricci tensors are equal to zero. 

One can always introduce a connection and another matrix valued function 
E so that L takes the form: 

L = -(<rV„V„ + £7) (8) 

Here V M is a sum of covariant Riemannian derivative with respect to metric g^ v 
and connection u^. One can, of course, express E and u> in terms of a M , b and 

w^^+f 1 ;). ( 9 ) 

E = b- gT{d^ + - u p Vp (10) 
For the future use we introduce also the field strength for u: 

Qfiu = d^LOy - d u u^ + [up, u v \ . (11) 

The connection will be used to construct covariant derivatives. The sub- 
script ;//... va will be used to denote repeated covariant derivatives with the 
connection u and the Christoffel connection on M. The subscript : a . . .be will 
denote repeated covariant derivatives containing uj and the Christoffel connection 
on the boundary. Difference between these two covariant derivatives is measured 
by the extrinsic curvature (JZJ). For example, E- a b = E- ao — L a bE- n . 

Let us define an integrated heat kernel for a hermitian operator L by the 
equation: 

K{Q,L,t):=Tr(Qexp(-tL))= d n Xyfgtr (Q(x)K(t; x; x; L)) , (12) 

Jm 
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where Q(x) is an hermitian matrix valued function, tr here is over matrix indices. 
For the boundary conditions we consider in this paper there exists an asymptotic 
expansion jH] as t — > 0: 

oo 

K(Q, L,t)~J2 L)t^ 2 . (13) 

According to the general theory jHj the coefficients aj.(Q,L) are locally com- 
putable. This means that each ak{Q,L) can be represented as a sum of volume 
and boundary integrals of local invariants constructed from Q, Q, E, the curva- 
ture tensor, and their derivatives. Boundary invariants may also include S, L ab 
and x- Total mass dimension of such invariants should be k for the volume terms 
and k — 1 for the boundary ones. 

At the moment several coefficients of the expansion ()13j) are known for the 
case of mixed boundary conditions © and matrix valued function Q (see [Zj for 
details of derivation; the formula (J4*§j) for a 4 was derived in [7j with additional 
restrictions L ab = and 5 = 0): 

a (Q,L) = (4n)- n / 2 [ d n x^tr(Q). (14) 

JM 

a 1 (Q,L) = ^(4 7 r)-("- 1 )/ 2 / d^xVhtrixQ)- (15) 

4 JdM 

a 2 {Q, L) = ^(4vr)- n / 2 d n x^/gti (6QE + QR) 

+ / d"" Wfctr (2QL aa + 12QS + 3 X Q-,n) ) ■ (16) 

JdM J 

a 3 (Q, L) = -J-(4 7 r)-("- 1 )/ 2 / d n - l xVhtr{Q(-24E + 2A X E X 

JdM 

+A8 X E + A8E X - l2 X :aX:a + ^X:aa ~ QX:aX:aX + 

+8 X R a i Kin 

+ 1925 2 + 96L aa S + (3 + 10 X )L aa L bb 
+(6 - A x )L ab L ab ) + Q ]n (9QS + 192S 2 ) + 24 X Q ;n „}. (17) 

For a scalar function Q and mixed boundary conditions the coefficients a 4 and 
as were already derived [TU] . 

2.2 (-function 

Zeta function of an operator L is defined by 

aoo = £-^, as) 

A 
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where the sum is over all eigenvalues of the operator L. The zeta function is 
related to the heat kernel by the transformation 

i r>+oo 

( L (s) = -— dtt 8 - l K{I,L,t). (19) 

r ( s ) Jo 

Residues at the poles of the zeta function are related to the coefficients of the 
heat kernel expansion: 

a k (I, L) = Res s=(n _ fe)/2 (r( S )a(s)) . (20) 

Here I is a unit matrix with a dimension of the matrix functions a M , b in Q.From 
([27Ijl it follows that 

a n (/,L) = a(0). (21) 

In Euclidean four dimensional space the zero temperature one-loop path in- 
tegral over the boson fields (ft = J2x Ca^a can be evaluated as follows (up to a 
normalization factor): 



(22) 



Here we introduced the constant fi with a dimension of mass in order to keep a 
proper dimension of the measure in the functional integral. (l(0) can be thought 
of as a number of eigenvalues of the operator L. For the operator L in the form 
the number of eigenvalues is infinite, so Cl(0) yields a regularized value for 
this number. 

The zero temperature one-loop effective action is defined then by 

W = -In Z = -~ IndetL + ^Cl(O) ln/i 2 = i&(0) + ^Cl(O) ln/i 2 = 

= ~G" 2s Cz(s))U (23) 

The term (l(0) In fi 2 = 0,4(1, L) ln/i 2 in the effective action W determines 
the one-loop beta function, this term describes renormalization of the one-loop 
logarithmic divergences appearing in the theory. 



3 Applications 

3.1 Free energy for boson fields 

A finite temperature field theory is defined in Euclidean space, since for boson 
fields one has to impose periodic boundary conditions on imaginary time co- 
ordinate (antiperiodic boundary conditions for fermion fields respectively). A 
partition function is defined by 

Z((3) =Tre- pH , (24) 
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where if is a hamiltonian of the problem and f3 = h/T. Let us choose the 
lagrangian density p in the form 

where r is an imaginary time coordinate and L is a three dimensional spatial 
part of the density in the form (|T|). The free energy of the system is defined by 

F(P) = --lnZ(P) = --ln(jV^ J D0exp(- J dr J d*x<j>p<j>)) , (26) 

the integration is over all periodic fields satisfying 0(r + (3) = <f)(r) (Np is a 
normalization coefficient). As a result the eigenfunctions of p have the form 
exjp(iruj n )<px, where u) n = 2imj (3 and L0 A = A0a- The free energy is thus equal 
to P] 

n=— oo A 

where we introduced ^-function 

+oo 



n=— oo A 



and the parameter /i with a mass dimensionality in order to make the argument 
of the logarithm dimensionless (also see a previous section). 
Then it is convenient to use the formula 

-1 /.+oo +°° 

coo = pn/ ^ M EE e> " +A ' ( 29 ) 

and separate n = and other terms in the sum. For n ^ terms we substitute 
the heat kernel expansion for the operator L at small t 

oo 

e~ M = K(I; L;t)^J2 <*k(I, L)t^ 2 (30) 

A k=0 

and perform t integration, then we arrive at the high temperature expansion 
(/3 — ► 0) for the free energy F: 



F/h = ~Cl{0) - ^Cl(O) Hp') + (4vr) 3 / 2 

Qs 1 1 f/3fi\ q 4 / , /fy 
/3 (4tt) 3 / 2 11 V 2vr / 16tt 2 V n 2tt 

n>5 v 



a 7r 2 oi (r(3) a 2 1 
'/F90 ~ /^47r 3 / 2 ~ /J 2 24 



(31) 
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Here = a,k(I,L), (r(s) = Ylti=i n ~ s ls a Riemann zeta function, Cl( s ) — 
\~ s is a zeta function of an operator L, 7 is the Euler constant. The first two 

terms on the r.h.s. of (j31j) follow from the n = term. 
The term 

-F— — T 4 (32) 



/? 4 90 ft 3 90" 

is the leading high temperature contribution to the free energy. 

The classical limit terms due to the equality Cl(0) = °3 can be rewritten as 
follows: 

T(-iG(0) + a<0)mJy=Tl>^. ,33) 

A 

The terms on the l.h.s. of ()33)1 yield a renormalized value of the terms on the 
r.h.s. of (|33jh since the sum on the righthandsight is generally divergent when 
the number of modes is infinite. 

The term with 04 determines the part of the free energy that appears due to 
one-loop logarithmic divergences and thus it depends on the dimensional param- 
eter as in the zero temperature case. 

3.2 Chiral anomaly in four dimensions 
for MIT bag boundary conditions 

Consider the Dirac operator on an n-dimensional Riemannian manifold 



(34) 



in external vector and axial vector fields. We suppose that and are 
anti-hermitian matrices in the space of some representation of the gauge group. 
crjf 7 ' is the spin-connection 1 . 

The Dirac operator transforms covariantly under infinitesimal local gauge 
transformations (the local gauge transformation is D — ► exp(— X)D exp(A)): 

5 X A^ = [An, X] 
hV^d^X+iV^X) 

D^D + [D,X] (35) 

and under infinitesimal local chiral transformations (the local chiral transforma- 
tion is D — > exp(iipy 5 )D exp(zy27 5 )): 

5 v An = d^ip + [V^cp], 
StpVp = -[A^ip], 

D^D + i{D,j 5 (p}. (36) 

1 The spin-connection must be included even on a flat manifold if the coordinates are not 
Cartesian. 
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The parameters A and (p are anti-hermitian matrices. 

We adopt the zeta-function regularization and write the one-loop effective 
action for Dirac fermions at zero temperature as 2 

W = -lndetS = -llndet D 2 = ^(0) + \ ln(// 2 )C S2 (0) , (37) 

where 

C 52 ( S ) = Tr(5- 2s ), (38) 

prime denotes differentiation with respect to s, and Tr is the functional trace. 
The following identity holds: 

Us) = TiA- s 6( A {s) = -sTr((6A)A- s - 1 ) . (39) 

Due to the identity (IBTfy 

8 x ( G2 (s) = - (2sTr([D, AID- 2 - 1 )) = -2s (Tr([D' 2s , X])) = 0, (40) 

so the effective action (pT7jl is gauge invariant, 5\W = 0. 

The chiral anomaly is by definition equal to the variation of W under an 
infinitesimal chiral transformation. Using (|3H|) we obtain: 

8Mb) = - (2i S Tr({A7V}S- 2s - x )) = -4is (Tr^V^)) , (41) 
and the anomaly reads 

A := 5 V W = ^Cg 2 (0) = -2Tr(z 7 V5- 2s )| s=0 . (42) 
The heat kernel is related to the zeta function by the Mellin transformation: 

POO 

Tr(i 7 Vir 2s ) = T(s)- 1 / dtf- 1 K(i'f<p,D 2 ,t). (43) 

Jo 

In particular, after the substitution of the heat kernel expansion (jlHjl into the 
formula (pfHj) we obtain 

A = -2a n (i 1 5 <p,D 2 ). (44) 
The same expression for the anomaly follows also from the Fujikawa approach 

u 

We impose local boundary conditions: 

n_v>U = o, n_ = i(i- 7 5 7n ) , (45) 

2 The one-loop effective action is proportional to Planck constant h, in what following we 
put h = 1 . 
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which are nothing else than a Euclidean version of the MIT bag boundary con- 
ditions jH]. For these boundary conditions III = and the normal component 
of the fermion current ip'^n^ vanishes on the boundary. 

Since D is a first order differential operator it was enough to fix the boundary 
conditions (f43j) on a half of the components. To proceed with a second order 
operator L = D 2 we need boundary conditions on the remaining components as 
well. They are defined by the consistency condition |12j : 

n_5^u = o, (46) 

which is equivalent to the Robin boundary condition 

(v n + s)n + ^u = o, n + = ^(i + 7 5 7n ) (47) 

with 

S = - l -Yi + L aa . (48) 

In the paper |7j the following expression for a coefficient a A {Q,L) with an 
hermitian matrix valued function Q and conditions (0), L a b = (flat boundaries), 
S = was obtained: 

a 4 {Q, L) = -i-(47r)- n / 2 { / d n xJgtr {Q(60£./ + GORE + 180£ 2 
360 I J M 

+30Q^O^ + + 5R 2 - IR^W + 2R tlupa R^)} 

+ / d n - 1 xVhtT{Q{30E. n + 30xE. n x + 90xE;n + ^0E ; nX 

JdM 

+18 X xAn + l2 X AnX + 18^ anX-a 

] + 54[ X :a, Han] + 30[*, ft a n:a] + l2R ;n + 30 X R-n} + 

+Q. n (-30E + 30 X E X + 90 X E + 90E X - 

-l8 X :aX:a + S0 X :aa ~ §X:aX:aX + 30 X i?) + 30 X Q;/ n } }• (49) 

To obtain the chiral anomaly in four dimensions 3 with MIT bag boundary 
conditions one has to calculate the coefficient a,4(Q,L) (|4*9"|) with L = D 2 , Q = 
irf<p and substitute it into (jSJ). We define V pu = d p V v — d v V p + [V^, V u ], A pv = 
D p A v - DyA^ D p A u = d p A v - T^ u A p + [V^AA, The anomaly contains two 
contributions: 

A = A v + A b . (50) 

3 In two dimensions (n — 2) the boundary part of the chiral anomaly with MIT bag boundary 
conditions is equal to zero [x]. 
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In the volume part 




-/ d l xJgti<p(-120[D ll V'»',A v \ 
Jm v 



+60 [D,A V , V» v ] - 60 D^D P D U A» + 120 {{D^A U , A"}, A 1 } 
+60 {D^, A U A U } + 120 A^D V A V A» + 30 [[A^, A u ],A pu ] 
+W<x {-45 % V pu V pa + 15 i A ph A pa - 30 z (V pu A p A a + AM"*"* 7 ) 
-120 % A p V up A a + 60 z A p A v A p A a } - 60 {D a A v )R va + 30 (D^A P )R 




only the DA — R terms seem to be new [Zj (for flat space it can be found e.g. in 



is new [Zj. It has been derived under the two restrictions: 5 = and L a b = 0. 
Note, that in the present context, the first condition (S = 0) actually follows 
from the second one (L ab = 0) due to (jlHj) . 
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